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Abstract 

We present a complete calculation of the Leading Order QCD corrections to the quark 
level decay amplitude for b —► S 77 and study their relevance for both the inclusive branch¬ 
ing ratio BR(B —> JQ 77 ) and for the exclusive decay channel B s —► 77 . In addition to the 
uncertainties in the short distance calculation, due to the choice of the renormalization 
scale, an appreciable uncertainty in both B s —> 77 and B —> Xg'j'y is introduced by the 
matrix element calculation. We also briefly discuss some long distance effects, especially 
those due to the rj c resonance for the inclusive rate. Finally, a brief analysis of the IR 
singularities of the two photon spectrum in the inclusive case is given. 


1 Introduction 

The radiative decays of the B meson are known to be very sensitive to strong interaction 
perturbative corrections as well as to the flavor structure of the electroweak interactions and to 
new physics beyond the Standard Model. In particular, both inclusive and exclusive processes 
induced by b — > sy have been studied in great detail [lj-|5j and two measurements already exist 
from the CLEO collaboration: BR(B —► A s y) = (2.32 ± 0.57 ± 0.35) x 10 ~ 4 and BR(B —■> 
K* 7 ) = (4.2 ± 0.8 ± 0.6) x 10“ 5 . 

Due to the impressive experimental effort which is being directed to the study of the physics 
of the B meson, we can be confident that much lower branching ratios will be measured in the 
future. Therefore it may be interesting to study processes induced at the quark level by a two 
photon radiative decay of the b quark, i.e. by b —> syy. 

The b —> S77 decay has received some attention in the literature [|6j-{|], because of the 
interest in the B s —> 77 exclusive mode. More recently, in Ref. 0 we focused on the study 
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of the inclusive B —> X ,.77 branching ratio. I 11 the pure electroweak theory, without QCD 
corrections but after the necessary kinematical cuts to isolate the contribution into hard photons 
are imposed, both branching ratios are found to be of order ICC'. There is at present an 
experimental upper bound on the BR(B S —> 77 ), namely BR(B S —> 77 ) < 1.48 x ICC 4 
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As we know from the study of b —> sy, the impact of QCD corrections on radiative B decays 
can be pretty dramatic. Therefore in this sequel, as we promised in ||, we now want to present 
the study of Leading Order QCD corrections to the quark level process b —► syy. We will use 
this result to predict the QCD corrected branching ratios for both the inclusive B —> AQ 77 
and the exclusive B s 77 mode. In both cases QCD corrections increase the branching ratio 
by 60% to more than 100%. On the other hand, the forward-backward asymmetry that was 
introduced in |Dj turns out to be very robust with respect to QCD corrections and always varies 
by less than 15%. 

I 11 order to motivate the interest of our perturbative calculation we will also comment about 
some relevant long distance contributions and devote particular attention to the effect of the 
r/ c resonance in the inclusive case. Moreover, we will see how some uncertainty for both the 
inclusive and the exclusive branching ratio is introduced at the level of the matrix element 
calculation, due to the dependence on m s . 

Finally, we will give in Appendix ^ the detailed description of the treatment of the IR 
singularities which arise in the spectrum of the two photons for B —»■ A s yy. 


2 Leading Order QCD corrections to b 


s 77. 


In this section we present the general structure of the leading Order QCD corrections to 
the quark level decay process b —> syy. We will give the expression for the amplitude 
A{b —> syy), including a complete resummation of the leading QCD corrections to all or¬ 
ders in (a s log(/x 2 /Mjy)) n . The result will be then specialized in the following sections to the 
calculation of the inclusive branching ratio BR(B —> AQ 77 ) and of the exclusive branching 
ratio for the decay B s —> 77 . 

We will discuss QCD corrections in the well established framework of electroweak effective 
hamiltonians with renormalization group improved resummation of QCD corrections. For a 
complete review of the subject see Ref. [|lTJ . The most general effective hamiltonian which 
describes radiative b —> s decays with up to three emitted gluons or photons is given by [|1^,|T3 
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( 1 ) 


Heff = 


v 1=1 


if*)°i 


where, as usual, Gf denotes the Fermi coupling constant and V l} indicates some Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element. In writing Eq. ([!]), we have used the unitarity 
of the CKM matrix and we have taken into account that for b —» s transitions V u bV* s -C 
V t bV t * ~ V c bV* s . The basis of local operators we use is obtained from the more general set of 
gauge invariant dimension five and six local operators with up to three external gauge bosons 


by appling the QED and QCD equations of motion [12|,[[3| and is expressed in terms of the 
following operators 


0 1 = 

(s a 7 /i Lc /3 )(c^7 / ,L6 Cf ) , 

o 2 = 

(s a ^Lc a )(cp'y t ,Lbp) , 

0 3 ,5 = 

{s a YLb a ) (qp r yn( L ,R)qp) , 

Ch,6 = 

(sa^Lbp) Y (9/37 fl (L,R)q a ) , 

q=u,...,b 

O7 = 

g 

1 Q^Sa^imhR + m s L)b a F^ , 

0 8 = 

l67T 2 Sa<J ^( m b R + m sL)tapbpG a lxv 


( 2 ) 


where the chiral structure is specified by the projectors L,R = (1 7 5 )/ 2 , while a and (3 are 

color indices. F IW and G“ u denotes the QED and QCD held strength tensors respectively, also 
e and g s stand for the electromagnetic and strong coupling constants. 

The Wilson coefficients C',;(/i) are process independent and their renormalization is deter¬ 
mined only by the basis of operators {Oj}. They depend on the renormalization scale /j which 
we will set eventually to fi ~ nib . This introduces an error in the theory that is quite significant 
when only Leading Order (LO) logarithms of the form (a s log(/x 2 /M^)) n are taken into account 
and gets appreciably reduced when also Next-to-Leading Order (NLO) logarithms of the form 
a s (a s log(/i 2 /M^))" are resummed. The LO result for the Wilson coefficients in Eq. ([[]) is now 
a well established result [[L]] and recently the authors of Ref. [Q provided us with the first NLO 
calculation. 

If we want to calculate the amplitude for b —> S 77 at LO we have to use the effective 
hamiltonian in Eq. ([lj) with LO Wilson coefficients and evaluate its matrix element for the 
b —> S77 decay at 0 (ct°). On the other hand, for a NLO result we have to use NLO Wilson 
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coefficients and include 0(a s ) corrections to the matrix element. 

In order to understand the impact of QCD corrections on this new class of rare radiative 
B-decays, we choose to perform our analysis including, for the time being, only LO corrections. 
Therefore we will take the LO regularization scheme independent Wilson coefficients from the 
literature 0 and will not consider explicitly the matrix elements due to the insertion of O 5 and 
Oq into the one photon and one gluon penguin diagrams. In fact these matrix elements are 
reabsorbed into the scheme independent definition of C 7 (p) and C$(fi) 

CV e » = CrM + Zl 8 • C(n) , (3) 

where C(/i) is the vector of C\(p ),..., C 6 (n), while the vectors Z 7 8 depend on the regulariza¬ 
tion scheme: they are zero in the ’t Hooft-Veltman (HV) scheme and non zero in the Naive 
Dimensional Reduction scheme (NDR) (see Ref. |3j] for details). I 11 our calculation, we use the 
Cf* effective coefficients, although we decide to drop the extra index to simplify the notation. 
We note that no new regularization scheme dependence enters in the calculation of the matrix 
elements for b —> S77 through the new class of penguin diagrams with two external photons. I 11 
fact, a finite scheme dependence in the matrix element can arise only as a result of the product 
of the UV pole part of a Feynman diagram (or set of diagrams) times some 0(e) evanescent 
Dirac structure of the diagram itself. However, as we will see, the new penguins with two 
external photons are UV finite at 0(a°). Therefore any difference between two regularization 
schemes can only give an unphysical 0(e) effect. We have performed the calculation of the 
following matrix elements in both the HV and NDR regularization schemes and, as expected, 
the results coincide. Therefore we do not specify any regularization scheme in the following 
discussion. 

The amplitude for the decay b{p) —> s(p') + 7 (^ 1 ) + 7 (^ 2 ) can be expressed as 

A = Ai = Ci(v)u s (p')Tru b (p)en(ki)6 u (k 2 ) , (4) 

*=1 v i=l 

where X t = V t bV t * and e^(ki) and e v [k 2 ) are the polarization vectors of the two photons. The 
Ci(p) coefficient are intended to be the LO ones, as explained before, while we have denoted by 
T[ IU the tensor structure of the transition amplitude induced by the operator Oj. The different 
T[ IU are obtained inserting the operators of Eq. (|2|) into the Feynman diagrams of Fig. [T]. 
according to the color and chiral structure of the operators themselves. In particular, one has 
to be careful when dealing with penguin-like operators 03 v ..,6 due to their more complicated 
flavor structure. The Tk u tensors can be summarized in a compact form as follows 
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Figure 1: Examples of Feynman diagrams which contribute to the matrix clement (s77|ZJ e ff \b)- 
The 1PI diagrams illustrate the two possible insertions of the operators Oi,...,Oq (double 
circled cross vertices), depending on their flavor, chiral and color structure, while the 1PR ones 
represent the insertion of 0 7 (cross vertices). Moreover q indicates a generic quark flavor. 


IT" = VQTUir , 

7T = QIkJVIT , 

Tg = { N c \_Q\ (Kd + K s + Kb) + Qu ( K U + Kcj] + Q\ (Kb + K s ) } , 

Tt = { \Qd ( K d + K s + K b ) + (fv u + Av c )] + N c Qd (. K s + Kfe)} Wjf , (5) 

rr = -N c [Ql(Kd + K s + K b )+Ql(K u + K c )\wr + Q 2 d[rn b W^ b R + m s W^ s L\ , 

Tr = -[Ql(K d + K s + K b ) + Ql(K u + K c )]Wr + N c Ql[m b W^R + m s W£ I ' s L\ , 

rr = Q d wr , 

where we note that there is no contribution from the chromo-magnetic operator Og at 0(a°). 
In Eq. (H) N c denotes the number of colors (N c — 3), Q u = 2/3 and Q d = — 1/3 are the up-type 
and down-type quark electric charges and m b and m s indicate the masses of the bottom and 
of the strange quark respectively. Moreover all the Tk u have been expressed in terms of only 
three tensor structures 

w.r = 

{ ki ■ k 2 

yy (Vi- Vi - <r (¥1 - m +l , (o 

Wr = T(7‘'¥»"¥i+U'*:27V + s'‘"Mi-^7"¥i-<¥fT'‘)(l-2K g ) + 

' mi 
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h u h^ 

4 («r - ^4-1«. 


<?) 


wr = 


2 [ 2p ■ k 2 
1 


h x ■ k 2 

k/iY(m b R + m s L)(J> - k/ 2 + m b )Y+ 


2p' ■ k 2 

and the analytic coefficients K q defined as 


Y(i> - k/x+ m a )^ff{m b R + m s L ) 


+ {{k\,p} <-► {k 2 ,v}) , 


Kq — - 

q 2 


Qo( Z q) 


1 1 1 dx 2 

= - H-/ - l 0 g(l - ZgZ + , 

2 Z q J 0 X 


(7) 


111 


for = 2fci • k 2 /m 2 q . In derivng Eq. (]5])-([7]) we have checked the analogous results given 
Refs. (T|,|T§ for the b —> syg decay and we confirm all of them. 

Finally, we observe that using the effective hamiltonian of Eq. ([[]) at // ~ ARv and in the 
absence of QCD corrections, we can reproduce the pure electroweak amplitude obtained in 
Refs. [;6j-|8|, as expected. Only two operators, 0 2 and 0 7 , contribute in this case. Their Wilson 
coeffcients at p = Mw read 


C 2 (M w ) — 1 , 

C 7 (M w ) = F 2 {x t ) - F 2 (x c ) , ( 8 ) 


where F 2 (xi ) is the Inami-Lim function for the on-shell bs'y vertex 0 


Q /"V»3 _ 0,y>2 _ _ PT/y.2 I ^'T* • 

F 2 ( Xi ) = 7 ^— 7 ^ log^ + ’ + * 


— l) 4 24(xj — l) 3 ' 

The corresponding matrix elements are given in Eq. (j^) and we can easily verify that 0 2 
reproduces the one particle irreducible part of the result of Refs. [6j-§ while O 7 is responsible 
for the one particle reducible part. 


3 Inclusive branching ratio for b — > S77. 

As already discussed in Ref. ||, the inclusive rate for B —> X s 77 can be described to a good 
degree of accuracy by the quark level process. We can therefore directly use the results of 
Section || to evaluate the square amplitude. For this purpose, we rewrite the amplitude as 


A = 


ie 2 G 


F , - 


y/2- 


7T 


A Mp') [F 2 Wr + F 5 (m b W£R + m s Wg a L) + F 7 W<T] u^e^k^h) , ( 10 ) 
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where the coefficients F t can be easily deduced from Eqs. (@)-@, and are 


F 2 = (N c Ci(fj) + C 2 (p))QlK c + 

C3 (aO {iV c \Q d (K d + k s + Kb) + Q 2 (k u + k c )] + Q d (K s + Kfi)} + 

C 4 (//) { [<5d(^d + K s + Kb) + Q 2 (k u + K c )] + N c Q 2 d ( K s + Kb)} — 

(N c C 5 (p) + Ce(fi)) \Q 2 d {K d + k s + Kb) + (k u + k c )] , (11) 

F 5 = (CM + N c CM)Ql , 

F 7 = C 7 {n)Q d . 

The square amplitude summed over spins and polarizations will then be given by 


\A\ 2 — - A t )' [| F 2 |"A 22 + !F 5 | 2 y4 55 + |F 7 | 2 74 7 7 + 2Re(F 7 F 2 *)A 2 7+ (12) 

2Re(F 5 F 2 *(l — 2K b ))A b 2b + 2Re(F 5 F 2 *(l — 2k s ))A 25 + 2Re(F 7 F^)A 57 ] , 

where the quantities A t] denote the contractions between the tensors W and W^ v . In order 
to give them explicitily we introduce the following notation 1 


s = 


2 ki ■ k 2 


t = 


2p ■ k 2 


u = 


2 p ■ k\ 


rn„ 

P=~2 


mi ml ml ml 


(13) 


which satisfy the relation: u + t — s = 1 — p. I 11 order to introduce a more compact notation, 
it can be useful to switch occasionaly to the (s, u, t) invariants, defined as s = s/(l — p), 
t — t/(l — p) and u — u/(l — p). In this framework the A i3 quantities are given by: 


A 22 

A55 

A b ’ s 

^25 

^27 

a§ 7 


2 [(1 — p) 2 — (1 + p) s ] , 

[16|k 6 | 2 + |(1 - 2 K b )s + AK b \ 2 + p (16|k s | 2 + |(1 - 2k s )s/p + 4k s | 2 )] (1 
16 Re (8 pK b K* s + s [K b - 2(1 + p)K b K * + />k*]} , 

±s(l-pTs) , 


(l + p)s + 


ps- 


+ 


ps 


(s — t)t (.s — u)u 



8(Kb + pK s )s - [4 p(n b + k s ) + s ((1 


2 K s ) + p(l - 2 Kb))] 


s + p) + 


(14) 


1 We decide to follow in our discussion the notation of Ref. 0 as closely as possible, which can be helpful 
for comparison. 
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A? 7 

with 


+ 


t(s — t ) u(s — u ) 

(1 + p) (1 — p)A 77 — 2Ajj 


, 4 ( 3 ) 

+ /1 77 5 


A 1 ) - 1 

/I77 — v 


A 2 ) - - 

/I77 — TV, 


A (3) _ 

/i 77 — 


2 w(u — 2 ) _ 2u — 1 -o 

1 + u + — p —^ + -P 

1 — u 1 — u 

t+ P _ xn p + (t 


+ (£ 


u ) 


1 + p 
1 — u 


{l-u} 


—2jz | (1 + p)(2 + wt) + — 1 — 

tu { 1 — p 


u) , 

2(1 + p) — tu 


(15) 


(1 -f)(l -11). 

We want to put particular emphasis on the structure of the A 77 part of the square amplitude 
because it will be a crucial ingredient in testing the cancellation of the IR divergences which 
appear in the calculation of the total rate. In fact, the total rate is obtained by integrating 


dr = 


5 D (p-p' -ki- k 2 )\A\ 


2 d <yD d^ D l ^k\ d^ 0 1) A)2 


(16) 


2m b (2n) 2D ~ 3 ^ ^ 2 p' 0 2u- v 2 uj 2 

over the physical phase space, where we have denoted by and u ; 2 the energies of the two 
photons. All the terms in \A\ 2 are both UV and IR finite except A 77 which gives origin to IR 
singularities upon integration over the phase space of the two photons. We chose to regularize 
the integrals working in D = 4 — 2e dimensions and to extract the existing IR singularities as 
poles in 1/e. These IR divergences originate when either uj\ — » 0 or oo 2 — 1 ► 0, and correspond to 
the well known IR singularities which arise in the bremstrahlung process when one or the other 
of the two photons becomes very softQ. I 11 this limit the b —> S77 decay cannot be distinguished 
from b —> sy and the two processes have to be considered together in order to obtain meaningful 
(i.e., finite) physical quantities. In fact, we have checked that the IR singularities which arise 
from the integration of A 77 over the phase space cancel exactly with the 0(a e ) virtual corrections 
to the b —> s'f amplitude (see Appendix |A|). Therefore (T(6 — > sy) + T(6 —> S77)) is free of IR 
singularities. 

This problem has already been studied in detail in order to take into account the 0(a s ) 
bremstrahlung corrections for b —> sy [|T4|, [lH] . However our point of view here is slightly 


different. In our case the bremstrahlung process is not considered as an 0(a e ) correction to 

2 We note that there are no collinear singularities so long as the mass of the external quarks are non zero. 

This gives origin to a non negligible dependence on m s and perhaps a more careful resummation of logarithms 
like log(m 2 ) in the rate should be implemented. We will discuss our concern with this problem later on. 


























the b —> sy amplitude, but as a different process: the decay of a b quark into an s quark 
plus two hard photons. Therefore, the endpoints of the spectrum of each photon (where the 
IR singularities are present) do not in fact correspond to the process of interest. In order to 
calculate the physical rate of interest we just have to impose a cut on the energy of each photon, 
which will naturally correspond to the experimental cut imposed on the minimun energy for 
detectable photons. 



0 0.2 0.4 0.6 0.8 1 


Xy 

Figure 2: The spectrum of the two photons including QCD corrections, normalized to the total 
QCD corrected rate for m s = 0.5 GeV. The two photons are defined as the photon of lower 
energy (solid) and the photon of higher energy (dashed). 

In Fig. we illustrate the spectrum of the two photons, defined as the photon of higher 
energy and the photon of lower energy. We obtain this spectrum requiring the energy of each 
photon to be larger than E 1 ™ 111 = 100 MeV and the angles between any two outgoing particles 
to be bigger than = 20°. This last constraint is not required analytically, but we think it 
is reasonable to exclude photons which are emitted too close to each other or to the outgoing 
s quark, in order to roughly incorporate the experimental requirements as we perceive them. 
Once the structure of the differential rate has been checked and the presence of IR singularities 
understood and treated, we can integrate Eq. ( [TUP numerically and study the impact of QCD 
corrections on the total rate as well as on different distributions. 

We find that QCD corrections enhance the rate by a factor of ~ 2 — 2.5, depending on the 
numerical parameters we use. In our evaluation we fix up, = 4.8 GeV, m c = 1.5 GeV, mt = 175 
GeV and |A*| = = 0.04. As far as m s is concerned, we use m s ~M K = 0.5 GeV. For this 
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set of parameters and fixing /x = m&, the branching ratio for B —> A s 77 goes from ~ 1.7 x 10~ 7 , 
without QCD corrections, to ~ 3.7 x 10 -7 when LO QCD corrections are included. We recall 
that we define the BR(B —> AT. 77 ) in terms of the scmileptonic branching ratio as follows 


BR(B —> AT 77 ) — 


r(6 —> S77) 


th 


x BR(B —> XJui 


exp 


(17) 
clvi ) at 


T(b —* clvi 

where no QCD corrections have to be included in the theoretical prediction of T(b 
this order in a s and we have used BR(B —> X c lv{) ~ 0.11 [|17] . 

In principle, rn s ~ Mx should be used in the phase space integration, while in the pertur¬ 
bative calculation of the amplitude one may need to replace it by the current mass m s ~0.15 
GeV. However this introduces spurious instabilities in the numerical Montecarlo integration 
over the phase space. Since the numerical results change little as we replace m s over the range 
0.15—0.5 GeV, we prefer to use the same value of m s in both cases. Thus, in order to simulate 
the physical phase space correctly we set everywhere. 

Moreover, we have to account for the scale dependence introduced by QCD corrections at 
the level of the Wilson coefficients. This makes a 25 — 30% uncertainty, as is the case for 
B —> AT 7 case. For the sake of completion, we also give the values of the Wilson coefficients 
we use in Table |I|. for three values of /i, fi = mb/2, m b and 2 mb respectively, and for m t = 175 
GeV and a s (Mw ) = 0.118. We will comment about further uncertainties introduced by long 
distance QCD effects in Section |5j. 


/J 

C\ 

c 2 

C3 

c 4 

C5 


CV 

m b j 2 

-0.324 

1.148 

0.015 

-0.033 

0.009 

-0.043 

-0.344 

rn h 

-0.234 

1.100 

0.010 

-0.024 

0.007 

-0.029 

-0.308 

2 mb 

-0.162 

1.065 

0.007 

-0.017 

0.005 

-0.019 

-0.277 


Table 1: Values of the regularization scheme independent LO Wilson coefficients for /i = mb/2, 
ji = mb, and ji= 2 mb, for mb = 4.8 GeV, m t = 175 GeV and a s (Mw) = 0 . 118 . 


I 11 order to better understand the dynamics of QCD corrections, let us classify the different 
contributions to the rate into one particle reducible (1PR) and one particle irreducible (1PI), 
as we did in Ref. || for the pure electroweak case. In the language of the effective hamiltonian 
of Eq. (f|) this corresponds to separating the cotribution of O 7 (which corresponds to the IPR 
part) from that of all the other operators. As we saw [jl|, the photon invariant mass distribution, 
dr /ds is dominated for low s by the IPR diagrams, while for larger s a non trivial m c -dependent 
contribution from the 1PI diagrams starts being relevant. The effect of QCD corrections is to 
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enhance even more the effect of 0 7 , as we conld expect from the dramatic effect that QCD 
corrections have in b —> s'y, while lowering the impact of the 1PI contribution, because of 
the new mixing with many different four-quark operators. In particular, the contribution of 
0 2 is suppressed by the destructive interference with 0\. We verified that the contributions 
of different operators to the angular distribution of the two photons are very similar to each 
other, also after QCD corrections have been included. 

On the other hand, as expected, the forward-backward asymmetry we introduced in pj 


Afb — 


r(cos 6> S7 > 0 ) - r(cos fl S7 < 0 ) 
r(cos60y > 0 ) + r(cos 6 Qy < 0 ) 
where 9 S1 is the angle between the s quark and the softer photon, is rather insensitive to 

QCD corrections, since the QCD corrections tend to cancel between the numerator and the 
denominator. In fact we find that QCD corrections affect %fb by no more than 15%, changing 
it from 0.71 (without QCD corrections) to 0.78 (with LO QCD corrections), despite the fact that 
the total rate changes by as much as 60% to 100%. Furthermore %fb is practically insensitive 
to the choice of scale in the LO Wilson coefficients, while the branching ratio varies as much 
as 30% with / 1 . On the other hand, this observable will clearly benefit from the enhancement 
induced by QCD at the rate level. Once the process is measured the possibility of measuring 
this new observable should give us another handle in testing our understanding of the theory 
and in differentiating the Standard Model from its extensions as already explained in ||. 


4 The exclusive decay B s 


77- 


Using the quark level amplitude in Eq. ([[(]) we can also estimate the rate for the B s 77 rare 
decay and evaluate the impact of QCD corrections on it. In order to calculate the matrix element 
of (LI);) for the B s —* 77 decay, we can work, for instance, in the weak binding approximation 
and assume that both the b and the s quarks are at rest in the B s meson. In the rest frame of 
the decaying B s meson we would have that 


h ■ k 2 = 


Ml 


Pb ■ ki = p b ■ k 2 = -m b M Ba 


Ps ■ k\ =p s -k 2 = -m s M Bs 


(19) 


where m b and m s must now be traeted as constituent masses. The problem can also be rephrased 
in the language of Heavy Quark Effective Theory (HQET), assuming that the velocity of the b 
quark coincides with the velocity of the B s meson up to a residual momentum of order Aqcd> 


i.e. p b = myv M + k To first approximation, the scalar products of Eq. (19) are replaced by 
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M I I 

ki-k 2 = —^- , p b ■ h = p b ■ k 2 = -m b M Bs , p s ■ k x = p s ■ k 2 = ~{M Bs - m b )M Bs ,(20) 

where we have used that p^ = — (p b — k\ — k- 2 Y. We can see that, to this order, Eqs. (pJ|) and (p0[) 
are compatible up to corrections of order (Aqcd / m b ) , if we assume rn s ps (M Bs — m b ) ~ Aqcd- 
Unless the HQET formalism is taken to beyond the leading order one cannot make a reliable 
distinction between the two predictions. Therefore, for concreteness, we give in the following 
the necessary matrix elements using the weak binding approximation. By further recalling that 


(0\s^ 5 b\B s (P Ba )) = -if Bs P } 


B s ) 


( 21 ) 


(0\s^ 5 b\B s ) = if Ba M, 


B s 


we obtain the following matrix elements for , W§\ u q and Wj V 


(0\WIT\ B s )e,(k 1 )e u (k 2 ) = 


1 

' 2 ‘ 

(OlWCIB.MfcMfc) = 

4 rn b m s 


/ 

1 — 2n b 8K b m b \ 


( 22 ) 


1 — 2k, 

+ I - b -iF, v F^ 

m b Mi J - 


1 — 2k s 8n s m s . „ 

+ -r^ ) 


m. 


Ml 


m b 
1 - 2k s 
m s 


iF 


(■ m b - m s ) 


F, P F^ - iF W F' U ' 


[m b + m s 

where f Bs denotes the B s meson decay constant. The amplitude A(B S —> 77 ) can therefore be 
expressed in terms of the only two tensor structures F IIU F IW and F IIV F IW 


A(B, 77 ) = —FF2\ t [a*F,„F^ + i A-F^) , 


(23) 


where F IW = 1/2 e lwprJ F pa . The coefficient A + and A of the CP-even and of the CP-odd term 
can be easely derived from Eq. fl22|) and read 


A + = 


A~ = 


4/b. 


4 /b - 


m b 

2 


-m Bs 

| (m 2 b - mg) 

rn h rn s 

—2F 2 — M Bs 


+ Mb 
M 2 B J ° Bs 


m x 


2 k s 8 K s m s 

+ ^IT 


F 7 

1 


2K b 1_ 


2k, 


m b 


rn x 


v ( m b + rn a ) 2 ' 

n 5 - 1 7 

m b m s 


(24) 

(25) 
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The QCD corrected coefficients F 2l F 5 and F 7 can be taken from Eq. (ED, while at 0(oP s ) 
they are simply given by F 2 = Q^k c C 2 (M w ), F 5 = 0 and F 7 = Q d C 7 (M w ), for C 2 (M W ) and 
C 7 {M w ) in d||). We notice that the terms proportional to F 7 in both A + and A~ are inversely 
proportional to m<|]. This is a clear signal of the relevance of non-perturbative effects to the 
evaluation of the matrix element for the decay rate of B s —> 77 . In the absence of a calculation 
of the matrix elements for this process which takes into account the higher order corrections 
in the HQET expansion, we can only give the perturbative prediction and try to estimate the 
theoretical error we have on that. Therefore we will use Eqs. ( |23| ) and (|24|) and vary m s in the 
range 300 <m s < 500 MeV. 

Let us first estimate the impact of QCD corrections on the rate 

M 3 4 * n o 2 

r ( B ^rr) = ^H^A,H|.47 + m 2 ) , (26) 

and on the ratio of the two coefficients A + and A~ 


R = 


l^f 

| A -| 2 


(27) 


As pointed out in Refs. ||p, the coefficients A + and A~ correspond respectively to photons 
with parallel (e(fci) • e[k 2 )) and perpendicular (e(fci) x e{k 2 )) polarization. The interest in the 
ratio R also crucially depends on the magnitude of the branching ratio itself and is therefore 
important to examine the impact of QCD corrections on both of t lieinQ 

I 11 the following we will use / Bs ~ 200 MeV, M Bs = 5.37 GeV, mb = 4.8 GeV, m c = 1.5 GeV, 
m t = 175 GeV and |A t | = = 0.04. Using the experimental life time of the B s meson, 

r s = 1.61 ■ 10 -12 s, we ffiid that the branching ratio BR(B S —» 77 ) goes, for m s = 0.5 GeV, from 
3.1 x 10 - ' without QCD corrections to 5.0 x 1CD 7 with LO QCD corrections, therefore increasing 
by about 62%. As far as A + and A~ are concerned, their ratio is substantially changed by the 
acion of QCD corrections. It goes from R = 0.28 without QCD corrections to R = 0.55 with 
LO QCD corrections. In fact at 0 (ck°) both A + and A~ depend on 1PR part of the amplitude 
(O 7 ) and only A~ is sensitive to the 1PI part (O 2 ). When we switch on QCD corrections, the 
contribution of 0 7 dominates and drives A + and A~ closer and closer. This effect is amplified 
by the cancellation which takes place in the 1PI sector, mainly among 0 2 and 0\. 

3 The matrix element of W^ v s does not scale as l/m s for small m s because also k s —> 1/2 for m s —► 0, therefore 
killing the l/m s terms. Moreover, the dependence on m s from this matrix element is very much suppressed by 
the smallness of the coefficient F 5 . 

4 One interesting implication of this is that A + /A~ can be used to construct a CP-violating observable 

which will pick up a dependence on Im(X t ) /Re(X t ) ~ OffiA 2 ), where r] and A correspond to the Wolfenstein 

parametrization of the CKM matrix. 
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The uncertainty on the perturbative calculation is dominated by the scale-dependence of 
the LO Wilson coefficients, which is around 25 — 30%. On the other hand, we estimate the 
uncertainty coming from non-perturbative QCD effects, i.e. from the calculation of the matrix 
element, to be of about 50%. Thus, attributing a 60% uncertainty to the central value (5x 1CT 7 ), 
we expect the branching ratio to be about 2—8 x 10 -7 . It would be very useful to have a more 
accurate calculation of these effects, perhaps by using HQET beyond the leading order, so that 
a more precise theoretical prediction can be obtained. Indeed it is not inconceivable that those 
corrections will further increase the branching ratio for B s —> 77 . 


5 Long distance QCD effects. 


As far as the B s —> 77 rare decay is concerned, as we discussed in the previous Section, we expect 
long-distance QCD corrections to be proportional to 1 /m s at the lowest order, introducing an 
uncertainty that asks for a more accurate computation of the matrix element. Other non 
perturbative effects could come from the formation of cc bound states in the decay process, i.e. 
from resonances. However in the B s —> 77 case these resonant states would be far off-shell and 
they are not likely to give a significant contribution to the rate (similar to the b —> sy case). 

The inclusive decay B —> W 77 is, in this respect, more problematic. In the region of 
invariant mass of the two photons around the rate is going to be dominated by the 

Tj c resonance, which subsequently decays into two photons, i.e. by b —> r] c s —*• 577. This could 
affect other regions of the spectrum and constitute a serious problem. Moreover, we remind that 
in the resonance region, the inclusive B —> A s yy decay cannot be approximated anymore by 
the quark level process, as is the case for B —> X s e + e~ [|J. In order to understand the relevance 
of our perturbative calculation we need to include the resonance at the amplitude level and to 
estimate how it affects the invariant mass ditribution, dT/ds, away from the resonance peak. 
This will allow us to select those regions of the spectrum which are free from major long-distance 
pollutions. In principle we should include in our analysis all the possible resonant channels. 
However the rj c resonance is dominant and is enough to provide us with an idea of the resonant 
effects. 

In order to model the contribution of the rj c resonance we need to provide an effective vertex 
both for the b — > s?/ c transition and for the r) c —> 77 decay that follows it. The bsr] c vertex can 
be derived from the amplitude for the b — > sr] c decay JT8|. Using the effective hamiltonian in 
Eq. (il|) and parametrizing the axial vector current matrix element 
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(0|c7 /i 75c|r/ c ) = -ifn c P% 

in terms of the decay constant f Vc ~300 MeV |0|, one gets0 


(28) 


(sr] c \H eS \b) = 


Cl + C 3 — C*5 + + C 4 — Cq) 


u s [—m s (l - 75 ) + m 6 (l + 75 )] u b 


x 


(29) 


For the values of the parameters used in this paper and taking the LO Wilson coefficients from 
Table [I], we can estimate BR(b —► ry c + anything) ~ 4 x 1CT 3 , more restrictive than the present 
experimental upper bound |jT7[| 


BR[b —> 7 ] c + anything) < 9 x 10 3 . (30) 

As far as the 7^77 vertex is concerned, we can assume the amplitude for r/ c —> 77 to be of the 
form 


A{r] c ^ 77) = iB~F , 
and use the experimental measurement 


(31) 


5 i?( 7 c -> 77) = 3 x 10 " 4 (32) 

to estimate \B~\ ~2.5—3xl0" 3 , for T Vc = 0.013 GeV and M nc ~3 GeV. The relative sign between 
the perturbative continuum and the resonant contribution can be determined via the same kind 
of unitarity arguments applied in Ref. jnij] to the b —> sip case. In fact, in the resonance region 
the perturbative amplitude is much smaller than the resonant one and therefore the relative 


sign between the two terms of the amplitude has to be positive, as in []T9 . 

The amplitude for the inclusive b —> S77 decay can now be written as the sum of a non¬ 
resonant, An r, and of a resonant, Ar part 


A(b -»• 577) = A nr + A r = A nr + u s (p')F R W^u b (p)e fi (k 1 )e^(k 2 ) , ( 33 ) 

where Anr, including LO QCD corrections, is given in Eq. (| 10 [) while Ar has been expressed 
in terms of the following coefficient and matrix element 
5 We assume simple factorization. 
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1 


Fr 

wg 


i—LB- 


4a P 


Ci + Cj, — C$ + —{C2 + C4 — Cq) 


S ~ M Vc + iT ^r M , 


r)c ±v± r)c 


2k^ l ' pa k lp k 2l7 [—m s (l - 75 ) + m b ( 1 + 75 )] 


(34) 


If we use Eq. (|I3|) to compute the invariant mass distribution of the two photons, we see that 
the effect of the resonance is very well localized around the resonance peak and does not affect 
in particular the region for s < 0.3. We can define in fact two regions, for 0.0 < s < 0.3 and 
for s > 0.5, in which the effect of the r] c resonance is practically negligible, as one can see in 

Fig. |. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

S 


Figure 3: The invariant mass distribution of the two photons in the presence of the rj c resonance, 
normalized to the total rate r to t = 5.7 x 10 -7 , as obtained for m s = 0.5 GeV. We show the pure 
non-resonant (solid), the pure resonant (dot-dashed) and the total distribution (dashed). The 
resonance peak is truncated in order to show the relevance of the different contributions both 
inside and outside the resonance region. 


Over these regions we can assume the validity of our perturbative calculation of Section 0 as 
well as of our previous studies of the various kinematical distributions for b —> S 77 decay [[]]. 
Disregarding in the perturbative calculation of Section 0 the contribution of the resonance 
region, which we conservatively define as 0.3 < s < 0.5, we find that the perturbative branching 
ratio is reduced by at most 14%. It would be very useful to verify experimentally that the effect 
of the 7 / c resonance in the B —> W .77 case is not so relevant, in comparison with what we know 
to be the case for B —> X s e + e~. In fact, if we consider the decay chain b —» s-0 followed by 
ijj —> e + e~ and use both experimental [[17]] and theoretical | 11 ,^(J inputs, we can estimate that 
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(35) 


r(6^sy.)r(V>^e+e-) _ 2 

r(& —> se + e~) 

while the analogous quantity for b —> srj c followed by rj c —> 77 amounts to 


h(> -»• s 7 1 c ) T { 7 1c -> 77) 


~ 6 . 


(36) 


r (6 —> S 77 ) 

This argument indirectly confirms the less dramatic impact that the 7 ] c resonance has on the 
invariant mass distribution of the two photons in the b syy decay. 
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Note added. While in the course of writing this manuscript, we became aware of the following 
two papers: G.Hiller and E.O. Iltan, |icp-ph/9704385| and C.-H. V. Chang, G.-L. Lin and Y.- 
P. Yao, |hcp-ph/9705345| , in which the problem of QCD corrections as it pertains only to the 
exclusive B s —> 77 decay is also discussed. We disagree with the Erst reference, where only the 
contribution of 0^ and O 2 is considered. We agree with all the results of the second reference, 
except for a few points that appear to be misprints. 


A Study of the IR divergences of the rate. 

In this appendix we want to show the explicit cancellation between the IR singularities arising 
respectively in the b —> syy rate from the bremstrahlung of a soft photon and in the 0 (a e ) 
virtual corrections to the b —> sy amplitude. This will confirm our understanding of the 
endpoints of the photon spectrum in the b —> syy decay. Our calculation is very similar to 
what can be found in Refs. [14, T3|] for the study of the gluon spectrum in b —> syg. Many 
results could be taken from there provided the different charge and color factors are adequately 
taken into account. We have indeed reproduced the calculation we report in this appendix and 
we can confirm^] a posteriori the results for b —»■ syg. 

Moreover, as we already explained in Sec. we are not going to include 0(a e ) virtual 
corrections to the b —> sy amplitude in the calculation of the rate for b —>■ syy. I 11 fact, we will 


6 In the course of these checks we came across a misprint in Eq. (34) of Ref. [jlj] . We are very grateful to the 
author for confirming this. The correct expression is given in Eq. (|3q). 
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just require the two photons to be hard, imposing a minimun energy cut. Therefore, in the 
present appendix we will consider only those aspects of the discussion which are necessary to 
show the cancellation of the IR poles. 

In the reaction b{jp) —> s(p') + 7 (^ 1 ) + 7 (^ 2 ), the spectrum of any of the two photons presents 
two sharp singularities in the vicinity of the endpoints, i.e. for x 7 —> 0 and for x 7 —> 1 , where 
we dehne x 1 = E 1 /E™ ax for E“ ax = (ml — m 2 s )/2/mi,. The variable x 7 corresponds in general to 
the reduced energy of a given photon. To make contact with the notation introduced in (fHf ), 
we can easily see that : 


x 7 i 


u 

1 ~ P 


u and x l2 


t 

1 ~ P 


= t 


(37) 


This singular behavior at the endpoints of the spectrum corresponds to the presence of IR 
singularities in the rate for b —> 577 , when the energy of one or the other of two photons goes 
to zero, i.e. when x y —> 0 (the energy of the photon under consideration) or x 7 —> 1 (the energy 
of the other photon). 

These IR singularities originate from the integration of the A 77 part of the square amplitude 
over the phase space of the two photons. As we can see from Eq. ( 0 ), A 77 is symmetric with 
respect to (u t), i.e. under the exchange of the two photons. Therefore the treatment of the 
two endpoints is symmetric. Given the spectrum of one photon, we will arbitrarily consider the 
endpoint x 7l —> 1. All our results will be valid in an analogous manner for the other endpoint, 
i.e. for x 72 —> 1 . 

Let us consider the contribution of 0^ only to the differential decay rate. Starting from 
Eq. (0) and working out the integration over the phase space in D — 4 — 2e dimension we get 


dV 7 
dt du 


= (1 - P) 


, 1 OL f 


2\2e 


r 0 F 7 


(1 _ 

2e 


4 7 T u 7 T(2 - 2e)[iu(l - f(i,u) 2 ) 1/2 } 


A 77 (t,ui) 


(38) 


where A 77 (t,u ) is given in Eqs. (|H|) and (|H]) and the function 


f(t,u) = 1 


2 (u + t — 1 ) 


(39) 


(1 — p)tu 

corresponds kinematically to the cosinus of the angle between the two photons in the rest frame 
of the b quark, when expressed in terms of the invariants u and t. Moreover we denote by To 
the quantity 

Ct'Ln IX.L _ 

(40) 


GW N 2 5 
T 0 = , —ml 


32vr 4 
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The origin of the singularity in u —> 1 becomes evident after we integrate over t and similarly 
for the singularity in t —*■ 1 when we integrate over ftf]. In particular, they are generated by the 


the first bracket in A$ (see Eq. (fTSp ), whose contribution, upon integration, reads 


where 


and 


,/pIR -i 

-M = (1 + p)(1 - p) ^r 0 F?c< 


2 (G^ + eG^) 

(1 - u) l+2t 


C e = 


(1 -p) (47T// 2 / ml) 

T (2 - 2 e)h 2e 


2\2e 


(41) 


(42) 


G {a) = 
G (b) = 


1 + 


P 


1 — u 


u + log(l - u) , 
1 - p 


P [2 + log(l — u)\u — 2 ^-— log(l — u) + ^ ^ 


1 — u 


(43) 


log 2 (l -u) - 2Li 2 (w) 


1 -p ov y 1 - p [2 

using the standard notation for the Spence function Li 2 (x) = — Jg 1 dt log(l — xt)/t. After the 
last integration over u, the IR singularity for u —> 1 appears as a pole in e, i.e. 


r£t-*(& -> S 77 ) = [ du d -^- = (1 + p)(l - p) 3 -—ToFy 1 
Jo 


du 


A -0 r 7 

4 7 T e 


2 + (j——log p 
1 ~ P 


+ 


, (44) 


where the dots indicate all other kinds of terms arising from the integration. An analogous 
singularity arises for £ —► 1 when we integrate the second term of A 2 7 (see Eq. (|T5|)) first over 
du and then over dt. Therefore the rate has a total IR singularity given by 


r®(t - f> 77 ) = 


dTy R ol OL e f, 1 

du— = (1 + p)(l - p) 3 -— T 0 E 7 2 - 


'0 


du 


2 7T 


2 + log P 

1 -p 


+ 0(1) , (45) 


where we have indicated with 0 ( 1 ) all the other non singular terms arising from the integration. 

We will now show that the same IR singularity, but with opposite sign, arises from the 
0 (a e ) virtual corrections to the b —> sy amplitude induced, of course, by the same operators 
O 7 . In this case, given the tree level vertex induced by O 7 , we have to consider both self¬ 
energy and vertex 0 (a e ) corrections in the renormalized theory, i.e. taking into account the 
wave function renormalization constants of the b and of the s quark. The choice of gauge for 


7 We could obtain this second singularity by looking, after we integrate over t, to the u —> 0 endpoint of the 
remaining integration. However, we prefer to use the symmetry between the two photons. 
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the photon is not relevant if the calculation is consistently performed (we checked the result in 
both the Feynman and the Landau gauge) and the final result reads 


r(“ e) (&-s7) 


(1 +P )(l- p) 3 r 0 F?(l + 2g)'‘ (1 - P)- 3 - 
r ( 7 0) + sr ( 7 ae) , (46) 


where 


r<. 0 | (f,^s 7 ) = (i + rt(i-p) 3 r 0 F| , 

= (i+p)(i-p) 3 roF|2/i t ,, r r ( ( 2 1 J 2 £ e ) ) (-LJii-p)-* ■ 


(47) 


K„_ = 


lOe 

4 7r 


(47t)T( 1 + e) 


— 2 + 


l+_p 

1 ~ P 


logp + 0(1) 


(48) 


It is now easy to verify that the pole terms cancel between Eq. (£F3) and Eq. (IM), such that 


T i 7 r (6 -»■ s 77 ) + 8T^\b -»■ s 7 ) = 0(1) . (49) 

In the previous discussion we may have disregarded terms of 0(e) when they do not happen to 
multiply a quantity containing 1/e poles and we have omitted all over a factor of (m^mj)/mj 2 
because it would not influence the cancellation of the IR poles. 
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